Precision tests of the Standard Model with leptonic and semileptonic 
kaon decays 

B. Sciascia ^ on behalf of the FlaviaNet Kaon Working Group * 
^Laboratori Nazionali di Frascati dell'INFN 

This paper presents the analysis of leptonic and semileptonic kaon decays data done by the FlaviaNet Kaon 
Working group, as described in [1]. Data include all recent results by BNL-E865, KLOE, KTeV, ISTRA+, and 
NA48. Experimental results are critically reviewed and combined, taking into account theoretical (both analytical 
and numerical) constraints on the semileptonic kaon form factors. We report on a very accurate determination of 
Vus as well as on many other tests of the SM which can be performed with leptonic and semileptonic K decays. 



1. Introduction 

In the Standard Model, SM, transition rates of 
semileptonic processes such as (f — > uHi/, with 
{u^) being a generic down (up) quark, can be com- 
puted with high accuracy in terms of the Fermi 
coupling Gp and the elements Vji of the Cabibbo- 
Kobayashi Maskawa (CKM) matrix. Measure- 
ments of the transition rates provide therefore 
precise determinations of the fundamental SM 
couplings. 

A detailed analysis of semileptonic decays of- 
fers also the possibility to set stringent constraints 
on new physics scenarios. While within the SM 
all cP — > u^iv transitions are ruled by the same 
CKM coupling Vji (satisfying the unitarity con- 
dition J2k l^ifeP — 1) and Gp is the same cou- 
pling appearing in the muon decay, this is not 
necessarily true beyond the SM. Setting bounds 
on the violations of CKM unitarity, violations of 
Icpton universality, and deviations from the V — A 
structure, allows us to put significant constraints 
on various new-physics scenarios (or eventually 
find evidences of new physics). 

In the case of leptonic and semileptonic K de- 
cays these tests are particularly significant given 
the large amount of data recently collected by sev- 
eral experiments: BNL-E865, KLOE, KTeV, IS- 
TRA+, and NA48. The analysis of these data 
provides precise determination of fundamental 

*WWW access at www.Inf.infn.it/wg/vus/; for sake of 
completeness and brevity for all references we refer to the 
Note written by the FlaviaNet Kaon WG [l]. 



SM couplings, sets stringent SM test almost free 
from hadronic uncertainties, and finally can dis- 
criminate between new physics scenarios. The 
high statistical precision of measurements and the 
detailed information on kinematical distributions 
have pushed a substantial progress on the the- 
ory side, in particular the theoretical error on 
hadronic form factors has been reduced at the 
1% level. 

The paper is organized as follows. First in 
Sec. [5] we present fits to world data on the lead- 
ing branching ratios and lifetimes, for K^^ Ks, 
and mesons. Sec. [3] summarizes the sta- 
tus of the knowledge of form factor slopes from 
Ki'i decays. The physics results obtained are 
described in Sec. [H in particular the measure- 
ment of \Vusf+{Q)\- Finally, to the special role of 
r(i^^=^2)/r(^;^2) ratio is devoted the Sec. O 

2. Experimental data: BRs and lifetime 

Numerous measurements of the principal kaon 
BRs, or of various ratios of these BRs, have been 
published recently. For the purposes of evaluating 
|Kis/+(0)|, these data can be used in a PDG-like 
fit to the BRs and lifetime, so all such measure- 
ments are interesting. A detailed description to 
the fit procedure and the references of all experi- 
mental input used can be found in Ref. [Ij. 

For Kl the results are given in table [1] while 
table [2] gives the results for K^. For the Ks^ the 
fit is dominated by the KLOE measurements of 
BR{Ks Tiev) and of BR{ti+ti-)/ BR{-k°tt°). 
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Parameter 


Value 


S 




0.4056(7) 


1.1 


BR(i^^3) 


0.2705(7) 


1.1 


BR(37rO) 


0.1951(9) 


1.2 


BR(7r+7r-7rO) 


0.1254(6) 


1.1 


BR(7r+7r-) 


1.997(7) X 10-3 


1.1 


BR(27rO) 


8.64(4) X lO-'^ 


1.3 


BR(77) 


5.47(4) X lO-'^ 


1.1 


TL 


51.17(20) ns 


1.1 



Table 1 

Results of fit to Kl BRs and lifetime. 



Parameter 


Value 


S 




63.57(11)% 


1.1 


BR(7r7rO) 


20.64(8)% 


1.1 


BR(7r7r7r) 


5.595(31)% 


1.0 


BR(Ke3) 


5.078(26)% 


1.2 




3.365(27)% 


1.7 


BR(7r7r07rO) 


1.750(26)% 


1.1 




12.384(19) ns 


1.7 



Table 2 

Results of fit to BRs and lifetime. 





Kl and K 


Measurements 


16 


xVndf 


54/13 (7 X 10-^) 


A'+ X 10^ 


24.9 ± 1.1 {S = 1.4) 


a:; X 103 


1.6 ±0.5 [S 1.3) 


Ao X 10^ 


13.4 ± 1.2 (5 = 1.9) 


p(a;,a';) 


-0.94 


p(A'+,Ao) 


+0.33 


p(A'|,Ao) 


-0.44 


/(^e°3) 


0.15457(29) 


^(^^3) 


0.15892(30) 




0.10212(31) 




0.10507(32) 




+0.63 



Table 3 

Averages of quadratic fit results for Kf.^ and K^^ 
slopes. 



present is better determined by measurements 
and by combining measurements and dispersion 
relations. Many approaches have been used, and 
all have been described in detail in [I]. Here we 
list only the averages of quadratic fit results for 
Kez and if^a slopes (Table [3]) used to determine 

|K.|/+(0). 



These, together with the constraint that the Ks 
BRs must add to unity, and the assumption of 
universal lepton couplings, completely determine 
the Ks leading BRs In particular, ^^.{Ks — > 
Tvev) = 7.046(91) x lO"''. For tks we use 
0.8958 X 10^1" s, where this is the non-CPT con- 
strained fit value from the PDG. 

3. Experimental data: Kg^ form factors 

The hadronic K n matrix clement of the vec- 
tor current is described by two form factors (FFs), 
f+{t) and fo{t). By construction, /o(0) = /+(0). 
In order to compute the phase space integrals 
we need experimental or theoretical inputs about 
the ^-dependence of FF. In principle, Chiral Per- 
turbation Theory (ChPT) and Lattice QCD are 
useful tools to set theoretical constraints. How- 
ever, in practice the t-dependence of the FFs at 



4. Physics results 

4.1. Determination of |T4s|/_|.(0) and 

\Vus\/\Vud\fK/U 

The value of |Kis|/+(0) has been determined 
from the decay rate of kaon semileptonic decays 
(see [1] for the detailed decomposition), using 
the world average values reported in previous sec- 
tions for lifetimes, branching ratios and phase 
space integrals. The results are shown in fig- 
ure [T] for Kl TTCv, Kl — > TTfiiy, Ks — > Trev, 
K^ — !■ Trei', K^ ir^v, and for the combi- 
nation. The average, |Ks|/-h(0) = 0.21664(48), 
has an uncertainty of about of 0.2%. The results 
from the five modes are in good agreement, the fit 
probability is 58%. In particular, comparing the 
values of |T4s|/+(0) obtained from K^^ and K^ 
we obtain a value of the SU(2) breaking correc- 
tion Sgjj,^) — 2.9(4)% in agreement with the 
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Figure 1. Display of |Kis|/+(0) for all channels. 



CHPT calculation (5fj^(2) = 2.36(22)%. More- 
over, recent analyzes on the so-called violations 
of Dashen's theorem in the kaon electromagnetic 
mass splitting point to <5fi/(2) values of about 3%. 

The test of Lepton Flavor Universality (LFU) 
between K^^ and iiT^a modes constraints a pos- 
sible anomalous lepton-flavor dependence in the 
leading weak vector current. It can therefore 
be compared to similar tests in r decays, but 
is different from the LFU tests in the helicity- 
suppressed modes and K12. The results on 
the parameter r^e = ^k'^s/k^J ^Kf^^/K^s is r^e = 
1.0043 ± 0.0052, in excellent agreement with lep- 
ton universality. With a precision of 0.5% the test 
in Ki3 decays has now reached the sensitivity of 
other determinations: Vf^eir) = 1.0005 ± 0.0041 
and r^e(7r) = 1.0042 ± 0.0033 2^ 

An independent determination of Vus is ob- 
tained from Ki2 decays. The most important 
mode is — > fj-'^v, which has been recently 
updated by KLOE reaching a relative uncer- 
tainty of about 0.3%. Hadronic uncertainties 
are minimized considering the ratio r(A'^ — > 
/i+i^)/r(7r+ l^^v). Using the world av- 
erage values of 'BK{K^ M**^) ^^'^ of 
given in Section [5] and the value of r(7r* 
^±v) = 38.408(7) ns-'^ from we obtain: 



IK. l/lKdl/x//. - 0.2760 ±0.0006. 

4.2. Theoretical estimates of /+(0) and 
Ik/ f-n 

The main obstacle in transforming these 
highly precise determinations of |y„s|/+(0) and 
\Vus\/\Vud\fKl f-K into a determination of |Ks| at 
the per-mil level are the theoretical uncertainties 
on the hadronic parameters /+(0) and fx/f-K- 
This hadronic quantity cannot be computed in 
perturbative QCD, but it is highly constrained 
by SU (3) and chiral symmetry. In the chiral 
limit and, more generally, in the SU{S) limit 
(m„ = rrid = vris) the conservation of the vec- 
tor current implies /+(0) = 1. Expanding around 
the chiral limit in powers of light quark masses 
we can write /+(0) = 1 -I- /2 + /4 + • ■ • where 
/2 and /4 are the NLO and NNLO corrections 
in ChPT. The AdemoUo-Gatto theorem implies 
that (/+(0) — 1) is at least of second order in the 
breaking of SU (3) This in turn implies that /2 
is free from the uncertainties of the 0{p^) coun- 
terterms in ChPT, and it can be computed with 
high accuracy: /2 = —0.023. The difhculties in 
estimating /+(0) begin with /4 or at 0{p^) in the 
chiral expansion. Several analytical approaches 
to determine /4 have been attempted over the 
years, essentially confirming the original estimate 
by Leutwyler and Roos. The benefit of these new 
results, obtained using more sophisticated tech- 
niques, lies in the fact that a better control over 
the systematic uncertainties of the calculation has 
been obtained. However, the size of the error is 
still around or above 1%, which is not compara- 
ble to the 0.2% accuracy which has been reached 
for |Ks|/+(0). 

Recent progress in lattice QCD gives us more 
optimism in the reduction of the error on /+(0) 
below the 1% level. Most of the currently avail- 
able lattice QCD results have been obtained with 
relatively heavy pions and the chiral extrapo- 
lation represents the dominant source of uncer- 
tainty. There is a general trend of lattice QCD 
results to be slightly lower than analytical ap- 
proaches. An important step in the reduction 
of the error associated to the chiral extrapola- 
tion has been recently made by the UKQCD-RBC 
collaboration. Their preliminary result /+(0) = 
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0.964(5) is obtained from the unquenched study 
with Np — 2 + 1 flavors, with an action that has 
good chiral properties on the lattice even at finite 
lattice spacing (domain- wall quarks). They also 
reached pions masses (> 330 MeV) much lighter 
than that used in previous studies of /+(0). The 
overall error is estimated to be 0.5%, which is 
very encouraging. 

In contrast to the semileptonic vector form fac- 
tor, the pseudoscalar decay constants are not pro- 
tected by the Ademollo-Gatto theorem and re- 
ceive corrections linear in the quark masses. Ex- 
panding fx I fn in power of quark masses, in anal- 
ogy to /+(0), /k/Ztt = I + r2 + ■■■ one finds 
that the C'(p'*) contribution r2 is already affected 
by local contributions and cannot be unambigu- 
ously predicted in ChPT. As a result, in the 
determination of JkIIt: lattice QCD has essen- 
tially no competition from purely analytical ap- 
proaches. The present overall accuracy is about 
1%. The novelty are the new lattice results with 
= 2-1-1 dynamical quarks and pions as light 
as 280 MeV, obtained by using the so-called stag- 
gered quarks. These analyzes cover a broad range 
of lattice spacings (i.e. a=0.06 and 0.15 fm) and 
is performed on sufficiently large physical vol- 
umes {m-^L > 5.0). It should be stressed, how- 
ever, that the sensitivity of fx / /tt to lighter pions 
is larger than in the computation of /+(0) and 
that chiral extrapolations are far more demand- 
ing in this case. In the following analysis we will 
use as reference value the MILC-HPQCD result 
Ik/U = 1.189(7). 

4.3. Test of CKM unitarity 

To determine | Kis | and | Vud \ we use the 
value |Ks|/+(0) = 0.2166(5) , the resuh 
\Vus\/\Vud\fK/U = 0.2760(6), /+(0) = 0.964(5), 
and fK/fyr — 1.189(7). From the above we 
find: IK^I = 0.2246 ± 0.0012 from Kes only, 
and \Vus\/\Vud\ = 0.2321 ± 0.0015 from Ki2 only. 
These determinations can be used in a fit together 
with the the recent evaluation of Vud from 0+ 
0+ nuclear beta decays: |Kd|=0.97418±0. 00026. 
This global fit gives Kd = 0.97417(26) and 
Vus = 0.2253(9), with xVndf = 0.65/1 (42%). 
This result does not make use of CKM unitar- 
ity. If the unitarity constraint is included, the 




Figure 2. Results of fits to \Vud\, \Vux\, and 

\Vus\/\Vud\- 



fit gives = 0.2255(7) and xVndf = 0.80/2 
(67%). Both results are illustrated in Fig. [2j 
The test of CKM unitarity can be also inter- 
preted as a test of universality of the lepton 
and quark gauge couplings. Using the results 
of the fit (without imposing unitarity) we ob- 

taim GcKM = + + IKbP]'^' = 

(1.1662±0.0004) X 10"^ GeV~^ in perfect agree- 
ment with the value obtained from the measure- 
ment of the muon lifetime: — (1.166371 ± 
0.000007) X 10"^ GeV"^ The current accuracy of 
the lepton-quark universality sets important con- 
straints on model building beyond the SM. For 
example, the presence of a Z' would affect the 
relation between Gckm and G^. In case of a Z' 
from S'O(IO) grand unification theories we obtain 
mz' > 700 GeV at 95% CL, to be compared with 
the niz' > 720 GeV bound set through the di- 
rect collider searches [2]. In a similar way, the 
unitarity constraint also provides useful bounds 
in various supersymmetry-breaking scenarios. 

4.4. sensitivity to new physics 

A particularly interesting test is the compari- 
son of the \Vus \ value extracted from the helicity- 
suppressed K12 decays with respect to the value 
extracted from the helicity-allowed Ke^ modes. 
To reduce theoretical uncertainties from fn and 
electromagnetic corrections in Ke2, we exploit the 
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ratio BR{Kf2)/BR{'Ki2) and we study the quan- 
tity 



Rl23 — 



0^ 



Within the SM, i?;23 = 1, while deviation from 1 
can be induced by non- vanishing scalar- or right- 
handed currents. Notice that in i?;23 the hadronic 
uncertainties enter through (/i<-//7r)//+(0). In 
the case of effect of scalar currents due to a 
charged Higgs, the unitarity relation between 
\Vud\ extracted from 0+ — > 0+ nuclear beta de- 
cays and I Vus I extracted from Kgs remains valid 
as soon as form factors are experimentally deter- 
mined. This constrain together with the exper- 
imental information of log C*^'^'^^'^ can be used 
in the global fit to improve the accuracy of the 
determination of i?i23, which in this scenario 
turns to be Rmll'^^i^, = 1.004 ± 0.007. Here 
(/ii://^)//+(0) has been fixed from lattice. This 
ratio is the key quantity to be improved in order 
to reduce present uncertainty on i?i23- This mea- 
surement of Ri23 can be used to set bounds on the 
charged Higgs mass and tan (3. Figure [3] shows 
the excluded region at 95% CL in the Mi/-tan/3 
plane. The measurement of BR(i? — > tv) can 
be also used to set a similar bound in the Mh- 
tan/3 plane. While B ^ tv can exclude quite 
an extensive region of this plane, there is an un- 
covered region in the exclusion corresponding to 
a destructive interference between the charged- 
Higgs and the SM amplitude. This region is fully 
covered by the K ^ fiv result. 

4.5. A test of lattice calculation 

The vector and scalar form factors f+.o{t) are 
analytic functions in the complex t-plane, except 
for a cut along the positive real axis, starting at 
the first physical threshold tth = {"m-K + itT'tt)'^, 
where they develop discontinuities. They are 
real for t < tth- Cauchy's theorem implies 
that f+.oit) can be written as a dispersive in- 
tegral along the physical cut where all possi- 
ble on-shell intermediate states contribute to its 
imaginary part. A number of subtractions is 
needed to make the integral convergent. Par- 
ticularly appealing is an improved dispersion re- 
lation recently proposed where two subtractions 




I 95% CL from K->|iv/7t->nv 
n 95% CL from B->tv 
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Figure 3. Excluded region in the charged Higgs 
mass-tan/? plane. The region excluded hy B ~^ 
TV is also indicated. 



are performed at t = (where by definition, 
/o(0) = 1) and at the so-called Callan-Treiman 
point tcT = {mj^ — m^). Since the Callan- 
Treiman relation fixes the value of scalar form 
factor at tcr to the ratio {fx/ f-rr)/ f+{^), the dis- 
persive parametrization for the scalar form factor 
allows to transform the available measurements 
of the scalar form factor into a precise informa- 
tion on {Jk/ fir)/ f+{Q), completely independent 
of the lattice estimates. Figure |4] shows the values 
for /-i-(0) determined from the scalar form factor 
slope measurements obtained using a dispersive 
parametrization and the Callan-Treiman relation, 
and /if/A = 1-189(7). from result on the FF 
slope using the dispersive parameterization The 
value of /+(0) = 0.964(5) from UKQCD/RBC is 
also shown. 



5. The special role of of V{Ke2) /y{K^2) 

The ratio Rk = r(X^2)/r(i^e2) can be pre- 
cisely calculated within the Standard Model. 
Neglecting radiative corrections, it is given by 



= 2.569 X 10"^ and re- 



n(0) ^ ^ 
K (m^-m2)2 

fleets the strong helicity suppression of the elec- 
tron channel. Radiative corrections have been 
computed with effective theories, yielding the fi- 



nal SM prediction i?^ — ^"-k 



?SM _ n(0) 



^^rad.corr. 
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Figure 4. Values for /+(0) determined from the 
scalar form factor slope using the Callan-Treiman 
relation and /k/U = 1.189(7). 



Figure 5. Exclusion limits at 95% CL on 
tan P and the charged Higgs mass Mff± from 
\Vus\Ke2/\Vus\Ki3 for different values of A 13. 



2.569 X 10-5 X (0.9622 ± 0.0004) = (2.477 ± 
0.001) X IQ-^. Because of the helicity suppres- 
sion within then SM, the Ke2 amplitude is a 
prominent candidate for possible sizable contri- 
butions from physics beyond the SM. Moreover, 
when normalizing to the K^2 rate, we obtain an 
extremely precise prediction of the i^e2 width 
within the SM. In order to be visible in the 
Ke2/ K^2 ratio, the new physics must violate lep- 
ton flavor universality. 

Recently it has been pointed out that in a 
supersymmetric framework sizable violations of 
lepton universality can be expected in K12 decays. 
At the tree level, lepton flavor violating terms are 
forbidden in the MSSM. However, these appear 
at the one-loop level, where an effective H^Ivt 
Yukawa interaction is generated. The non-SM 
contribution to Rk can be written as i?^^^ « 



dSM 



1 



M 



H± 



131 



' tan^ /3 



where 



A 13, the lepton flavor violating coupling , being 
generated at the loop level, could reach values of 
0{1Q~^). For moderately large tan/3 values, this 
contribution may therefore enhance Rk by up to 
a few percent. 

Experimental knowledge of Ke2 / Kf^2 has been 
poor so far. The current world average of Rk = 



BR(is:e2)/BR(if^2) = (2.45 ± 0.11) x 10"^ dates 
back to three experiments of the 1970s [2] and 
has a precision of about 5%. Three new pre- 
liminary measurements were reported by NA48/2 
and KLOE (see [T] for details). Both, the KLOE 
and the NA48/2 measurements are inclusive with 
respect to final state radiation contribution due 
to bremsstrahlung. Combining these new results 
with the current PDG value yields a current world 
average of Rk = (2.457 ± 0.032) x 10-^ with a 
relative error of 1.3%, a factor three more precise 
than the previous world average. This value is 
in very good agreement with the SM expectation 
and gives strong constraints for tan (3 and Mfj± , 
as shown in Fig. [51 For values of A13 w 5 x 10"^ 
and tan/3 > 50 the charged Higgs masses is 
pushed above 1000 GeV/c^ at 95% CL. 
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